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Abstract 

In this paper we study tiie finite temperature SU(2) gauge tiieory 
with staggered fermions for non-zero imaginary and real chemical po- 
tential. The method of analytical continuation of Monte Carlo results 
from imaginary to real chemical potential is tested by comparison with 
simulations performed directly for real chemical potential. We discuss 
the applicability of the method in the different regions of the phase 
diagram in the temperature - imaginary chemical potential plane. 



1 Introduction 



Understanding the QCD phase diagram in the temperature - chemical po- 
tential (T, fi) plane has many important implications in cosmology, in as- 
trophysics and in the phenomenology of heavy-ion collisions (for a review, 
see PJ). In particular, it is extremely important to localize with high accu- 
racy the critical lines in this plane and to determine if the transition across 
them is first order, second order or just a crossover 

The formulation of the theory on a space-time lattice offers a unique tool 
for such task. However, at non-zero chemical potential the determinant of the 
fermion matrix becomes complex, thus preventing to perform standard Monte 
Carlo simulations. One way out is to take advantage of physical fluctuations 
in the thermal ensemble generated at /i = in order to extract from it 
information at (small) non-zero fi, after suitable reweighting HI El- 
However, it is not possible to know a priori to what extent the reweighted 
ensemble overlaps the true one at non-zero fi} 

An alternative approach is based on the idea to perform numerical simula- 
tions at imaginary chemical potential jH] , for which the fermion determinant 
is again real, and to analytically continue the results to real fj, P HDl HH 1121 
IT^.^ In practice, the method consists in fitting the numerical data obtained 
for imaginary chemical potential n = i^i by a polynomial, which is then 
prolongated to real fj,. The advantage of this method lies in the fact that its 
applicability is not restricted by large lattice volumes, since it makes no use of 
reweighting. Therefore, size scaling analysis can be made after the analytical 
continuation to clarify the nature of the transition across the critical lines. 
On the other side, a limitation comes from the presence of non-analyticities 
in /i, thus implying, as we will see later on, that the method works practically 
only for yU < T. 

In particular, if non-analyticities are singled out for imaginary fi and even 
if numerical data are obtained inside a region of analyticity in the (T, fij) 
plane, it is necessary that during the "rotation" from /i = ifij to real fi = fiR 
no other non-analyticities are met. It is understood that non-analyticities 
show up only in the infinite volume limit and that all observables are analytic 

^In the (2+l)-dimcnsional U(l) Gross-Neveu model a clear disagreement was evi- 
denced p in the /i-dependence of the chiral condensate and of the number density be- 
tween the determinations from reweighting and from standard Monte Carlo techniques. 
Of course, this does not imply the ineffectiveness of the method in other theories and in 
different temperature regimes. 

^Some interesting arguments about the QCD phase diagram for imaginary chemical 
potential have been discussed in Ref. |14j . 
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on finite lattices. However, the presence of non-analyticities in tlie tliermo- 
dynamical limit is indicated on (large) finite volumes by sharp changes in the 
behavior of physical observables, which is practically difficult to interpolate 
by polynomials, if not by including a large number of terms. The papers 
where the method of analytical continuation has been used so far limited 
the attention to regions of temperature and chemical potential where this 
situation does not occur. 

Of course, the compatibility between the results from two very different 
approaches, such as reweighting and analytical continuation, would give rea- 
sonable confidence on the validity of both. Such consistency has been shown 
in the case of 4-fiavor QCD in Ref. [TT] . 

Nevertheless, a test of both methods by comparison with direct determi- 
nations from standard Monte Carlo is interesting, especially because it could 
shed light also on the applicability outside the regions usually considered in 
literature. For this purpose, it is necessary to turn to a theory where direct 
Monte Carlo simulations in presence of real chemical potential are feasible. 

It is well known that the determinant of the fermion matrix of 2-color 
QCD or SU(2) gauge theory is real even in presence of a real chemical po- 
tential (see, for instance [S]). Moreover, although being very different in the 
structure of the spectrum jT^, 2-color QCD shares many similarities with 
true QCD, therefore represents the best candidate for the mentioned test. 

In this paper we adopt 2-color QCD to compare determinations of the 
Polyakov loop and of the chiral condensate from the method of analytical 
continuation with results from direct Monte Carlo simulations. In particular, 
we characterize different regions of the phase diagram in the (T, fij) plane 
according to the applicability of this method. Our findings can be helpful to 
apply and to understand the results of the method of analytical continuation 
in the theory of physical interest, i.e. in true QCD. 

The paper is organized as follows: in Section 2 we describe the structure 
of the phase diagram of 2-color QCD in the (T, fij) plane and discuss the 
implications on the method of analytical continuation; in Section 3 we present 
the numerical results; in Section 4 we draw some conclusions and anticipate 
some future developments. 

2 The phase diagram of 2-color QCD 

In order to draw the phase diagram of 2-color QCD we make use of known the- 
oretical and numerical results on the behavior with temperature and chemical 
potential of the Polyakov loop L and of the chiral condensate (ipip). These ob- 
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Figure 1: Schematization of the phase diagram in the (T, 6) plane according 
to Roberge and Weiss [TB] . 

jects represent true order parameters only in two limiting cases: the Polyakov 
loop is the order parameter of the confinement- deconfinement transition in 
the limit of infinite quark masses, while the chiral condensate is that of the 
chiral transition in the limit of vanishing quark masses. However, there is a 
broad numerical evidence that they both exhibit a rapid variation in corre- 
spondence of certain values of the parameters of the theory, which therefore 
are taken as critical ones. 

We start the discussion from the case of imaginary chemical potential, fi = 
ifii. Roberge and Weiss have shown [TH] that the partition function of any 
SU(N) theory is periodic in the parameter 6 = fii/T (the Boltzmann constant 
is set equal to one) as consequence of a remnant of the Z(N) symmetry 
valid in absence of fermions. Moreover, they have shown by weak coupling 
perturbation theory that above a certain temperature Te the free energy and 
the Polyakov loop show cusps or discontinuities for 6 = 27r(A; + 1/2) /N, with 
k = 0, 1, . . .. This suggests that there are first order vertical critical lines 
in the {T,6) plane located at = 27r(A; + 1/2) /N and extending from Te 
to infinity (see Fig. 1). They argue also that below Te, there is no phase 
transition for any value of 6, this being true also in the physical case of 
e = fi = 0. 

The latter conclusion cannot be true for all values of Uf and of the fermion 
masses. In the chiral limit of QCD, there are arguments according to 
which a chiral phase transition for /i = does exist and is first order for 
Uf > 3 and second order for Uf = 2. For SU(2), where a strong coupling 
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Figure 2: Schematization of the phase diagram in the (T, 6) plane after the 
inclusion of the chiral critical lines. 

analysis could be applied, according to Ref. [TH] there is a second order chiral 
phase transition for Uf = 1 and nf = 2. For non- vanishing masses, this phase 
transition probably becomes a crossover. Furthermore, a numerical study on 
a 12^ X 4 lattice in SU(2) with nf = 8 and am = 0.07 has shown a peak in the 
susceptibility of the chiral condensate at Pc — 1.41, suggesting the possibility 
of a chiral phase transition near this point |19j . 

In view of all this, the worst case for the purposes of the analytical con- 
tinuation (but the most interesting from the point of view of the present 
study) is that the Roberge- Weiss (RW) critical line at 9 = n/N could be 
prolongated downwards and bent up to reaching the ^ = axis in corre- 
spondence of a critical temperature Tc (see also [TH [I2|)- For the reasons 
mentioned above, this prolongation of the RW critical line will be referred 
to in the following as the chiral critical line |^ . The RW periodicity 
and the symmetry for /i — — /i, resulting from CP invariance, imply that 
knowing the expectation value of an observable for T and 6 inside the strip 
{0 < < tt/N, < T < cxd} in the {T,6) plane is enough to fix it on 
the remaining part of the plane. Applying this argument also to the critical 
lines, we end up with the qualitative structure of the phase diagram drawn in 
Fig. 2. The actual nature of the transition across the branches of the critical 
lines of Fig. 2 which lie below Te is to a large extent unknown. Other sce- 
narios of phase diagrams are possible, depending on and on the fermion 
masses, such as, for instance, that where a chiral critical line emanating from 
Tc ends somewhere before joining the vertical lines |13j . 
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Inside the mentioned strip we can single out two regions where the (non- 
perturbative) determinations of an observable can be interpolated by a Taylor 
series in 6'^ = [fij/T)'^: the region above and the region below the chiral 
critical line. This Taylor series can be analytically continued to real n = 

by the "rotation" fii — > —ifiR, if no non-analyticities are met during 
this rotation. Such a condition is not a priori guaranteed. It is evident, 
that the crucial point is how the critical line in the (T, 6) plane, implicitely 
defined by T = f{fij/T), evolves during the rotation from imaginary to 
real /x. A reasonable assumption is that, at the end of the rotation, the 
critical line in the plane with real fi is given by T = f{—ifiji/T) = g{^ji/T), 
with g a monotonically decreasing function in ^r/T. Another assumption, 
which needs confirmation by numerical tests, is that the critical line evolves 
smoothly during the rotation, i.e. it spans in the {T, ^ir/T, fii/T) space a 
2-dimensional surface, connecting the critical line in the (T, fii/T) plane with 
the one in the {T,fiR/T) plane, which does not possess bumps or ridges or 
valleys. 

If these assumptions are true and the scenario of Fig. 2 holds, we can iden- 
tify three regions in the strip {0<^^<7r/A^, 0<T<oo}of the (T, /i//T) 
plane according to the possibility to apply the analytical continuation: 

• T > Te'- analytical continuation possible for < fii/T < 7c/N; 

• Tc < T < Te'. analytical continuation possible for < Hi/T < fi*j/T, 
where fi*j/T is the solution of T = f{fi*j/T), i.e. the value of fij/T on 
the chiral critical line at the given temperature; 

• T < Tc'. analytical continuation possible for < fii/T < —ifi*^/T, 
where ^i*r/T is the solution of T = 5f(/i^/T), i.e. the value of ^ir/T on 
the critical line of the {T,fiR/T) plane at the given temperature. 

So far the method of analytical continuation has been applied, implicitely 
making the assumptions we have mentioned, to the chiral condensate in the 
region T > Te for QCD with rif = 4 JT] and to the critical line itself for 
QCD with Uf = 4 [TI] and with nj = 2, 3 ^HI- In the case of QCD with 
Uf = 4, the analytically continued critical line has been found in agreement 
with the determination of Ref. 

In the next Section we will compare the scenario sketched in the above 
list with numerical simulations in all the three temperature regions. As a 
test-field theory we use 2-color QCD with nf = 8 and as observables we 
adopt the Polyakov loop L and the chiral condensate (ipip)- 
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3 Numerical results 



We consider the formulation of 2-color QCD or SU(2) gauge theory on a x 
A'^ lattice with staggered fermions for non-zero temperature and chemical 
potential of the baryonic density. 

The finite temperature is realized by compactifying the time direction and 
by imposing (anti)periodic boundary conditions on (fermion) boson fields. 
The connection between the temperature T and the lattice size in the time 
direction is T = l/{Nra), where a = a(/3) is the lattice spacing and P = A/g"^ 
is the bare coupling constant. The chemical potential is introduced by the 
replacements [20] 

f/t ^ e-'"'Ul, (1) 

where Ur and U}. are the gauge link variables in the forward and backward 
time direction, respectively. When /x is purely imaginary these replacements 
amount to add a constant background U(l) field to the original theory. For 
SU(2) the fermionic determinant, appearing in the partition function after 
integration of the fermion fields, is real for any complex value of the chemical 
potential fi, thus allowing the Monte Carlo importance sampling without any 
limitation. 

We adopted in this work the same hybrid simulation algorithm of Ref . PH] • 
Since for /i 7^ it is not possible to use the even-odd partitioning, we were 
forced to simulate eight degenerate continuum flavors. We performed all the 
simulations on a 8^ x 4 lattice, setting the quark bare mass to am = 0.07 
and the microcanonical time step to dt = 0.02. We made refreshments of the 
gauge fields every 5 steps of the molecular dynamics and of the pseudofermion 
fields every 3 steps. In order to reduce autocorrelation effects, we made 
"measurements" every 50 steps and analyzed data by the jackknife method 
combined with binning. 

Before presenting the numerical results, it is convenient to redraw the 
phase diagram discussed in the previous Section in units more suitable for the 
lattice. In particular, on the vertical axis we put the inverse coupling (3 which 
is connected to the temperature T by a monotonic function in the region of 
interest (increasing (3 means increasing temperature). On the horizontal axis 
we put the chemical potential in lattice units, fi = afi. The qualitative 
structure of the phase diagrams in the fij) and in the fiji) planes 
is the same as in the (T,^j/T) and in the {T,fiji/T) planes. In the new 
units, however, the RW transition lines in the (/?, fij) plane appear at fii = 
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2Ti{k + l/2)/{NNr). For = 2 the first RW transition line is located at 
fij = 7i/{2Nt-). Of course, the scenario illustrated at the end of the previous 
Section must be reformulated in terms of P, Pe and Pc- 

An estimate of the value of Pc can be deduced from the hterature. Indeed, 
in Ref. ^H], a study of the chiral condensate on a 12^ x 4 lattice for bare quark 
mass am = 0.07 lead to the result Pc = 1.41 ±0.03. Since we have a different 
lattice, we expect a small difference between our Pc and the determination 
of Ref. [in]. However, for our purposes we can safely adopt this value of Pc- 

As a first step we took two P values, one (much) above Pc and likely above 
Pe, the other below Pc in order to possibly single out a different behavior in 
fii of lattice observables across /t/ = 7r/(2A^^). We took /?=1.80 for the higher 
P value and P = 1.35 for the lower and determined the Polyakov loop and the 
chiral condensate for fij ranging between and the value corresponding to 
the second RW critical line, i.e. 37r/(2A'^^). The results are shown in Figs. 3 
and 4. We can clearly see the RW periodicity in all cases; however, for 
P = 1.80 we have discontinuity for the Polyakov loop and cusp for the chiral 
condensate across the first RW critical line, while for P = 1.35 the behavior 
is smooth for both observables. We have verified that the overshooting of 
the first RW critical line seen in Fig. 3 is due to an hysteresis effect, which 
confirms the first order nature of the transition. 

In order to proceed with the test of the method of analytical continuation 
we need also an estimation of Pe- There is a nice procedure to do it, suggested 
in Ref. It relies on performing simulations at a critical value of fij, say 
fij = 7r/(2A^^), and on changing P from a very high value downwards. If 
the zero field configuration is chosen as the starting one, the system remains 
always in the phase on the left of the RW line until P > Pe, if the flip to 
the other phase is prevented by a large space volume. Due to our limited 
computational power, we are constrained to a space volume which is too 
small to prevent flips between the two phases. An alternative procedure to 
determine Pe could be to study the behavior of the Polyakov loop (or of its 
phase) and of the chiral condensate across the RW critical lines for different 
P values and to see when the discontinuities or the cusps are smoothed out. 
Another possibility could be to exploit the hysteresis effect in the Polyakov 
loop or in the chiral condensate determinations across the RW critical line. 
The procedure consists in building two sequences of expectation values of the 
Polyakov loop; in the first (second) sequence the starting configuration of any 
new run is taken to be a thermal equilibrium configuration being always in 
the phase on the left (right) side of the first RW transition. Across the RW 
critical line the comparison of these two sequences of data shows a hysteresis 
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effect as long as /9 > (3e- This effect disappears below f3E- Of course, both 
the above procedures demand for a high statistics and, so far, we can only 
say that j3E is bounded between about 1.53 and about 1.57. This level of 
accuracy is, however, enough for the purposes of the present study. 

We turn now to the main part of this work, i.e. the test if the scenario 
presented at the end of the previous Section holds. We selected three values 
of (3 (1.90, 1.45, 0.90), each representative of one of the regions (3 > /3e: 
Pc < P < Pe and /3 < /3c- For each /3 we have performed simulations for 
both real and imaginary chemical potential, varying /t/ and jiR between 
and 7r/(2A^T-) and have determined expectation values of the Polyakov loop 
and of the chiral condensate. Then, we have interpolated the data obtained 
for imaginary /i with a truncated Taylor series of the form a + b fij + c ftj. 
After "rotating" this polynomial to real chemical potential, thus leading to 
a—b fijt+c we have compared it with the determinations obtained directly 
for real /i — /ir. 

For P — 1.90 we have interpolated with the fourth order polynomial 
the data obtained for < /t/ < 7r/(2A^^) and have found that the rotated 
polynomial perfectly interpolates data obtained directly for real fj, — /ir in 
the same range (see Fig. 5). The uncertainty in the rotated polynomial has 
been obtained by propagating in standard way the uncertainty on the fitted 
coefficients a, b and c. We have verified that the same situation occurs for 
other values of (3 larger than Pe- From this outcome we can conclude that 
the analytical continuation works in the region P > Pe for all the /i/'s before 
the first RW critical line. 

For P = 1.45 we are in the region P^ < P < Pe for which, if the scenario 
presented at the end of the previous Section holds, we expect that by varying 
fij we meet the chiral critical line, while no transitions should be met at the 
same P by varying fiR. Therefore this time we have interpolated with the 
fourth order polynomial the data obtained for the real chemical potential 
< /ii? ^ vr/ (2A^^) and have rotated this polynomial to its counterpart in /i/. 
The comparison with data obtained directly for imaginary chemical potential 
shows agreement for fii below ~ 0.18, which can be taken as an estimate of 
the chiral critical value in the {P,fii) plane at the given P (see Fig. 6). 

For P — 0.90 < Pc, we expect analyticity for all possible values of /i/, 
while the chiral critical line in the (/3, P^r) plane is met by varying p,R at 
the given value of p. Therefore we have interpolated with the fourth order 
polynomial the data obtained for < fij < 7r/{2Nr) and have rotated this 
polynomial to its counterpart in ftR. The comparison with data obtained 
directly for real chemical potential shows agreement for jjR below ~ 0.12, 
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which can be taken as an estimate of the critical value in the (/?, jln) plane 
at the given P (see Fig. 7). 

We can conclude from these results that the scenario of the end of the 
previous Section is consistent, thus giving support to the assumptions un- 
derlying it. 



4 Conclusions and outlook 

In this paper we have presented a possible structure of the phase diagram 
in the temperature - imaginary chemical potential in 2-color QCD and have 
argued how this structure can be continued to real chemical potential. Then, 
we have shown that the scenario we have described is consistent with Monte 
Carlo numerical determinations of the Polyakov loop and of the chiral con- 
densate for real and imaginary chemical potential. 

In particular, we have found that the method of analytical continua- 
tion, largely exploited now in the literature for the physically interesting 
case of QCD, works fine within the restrictions posed by the presence of 
non- analy t icit ies . 

High precision determinations were not the aim of this paper. We limited 
ourselves to give convincing numerical support to the proposed scenario, 
postponing massive calculations to forthcoming works. 

We plan, moreover, to further develop the present ideas by (a) determin- 
ing the location of the critical lines for real and imaginary chemical potential 
by the standard method of susceptibilities, (b) verifying that the critical lines 
for real chemical potential can be obtained from analytical continuation of 
those for imaginary chemical potential, (c) following the "evolution" of the 
critical lines during the rotation from imaginary to real chemical potential 
by simulations with complex chemical potential. 
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Figure 3: Polyakov loop (above) and chiral condensate (below) vs fij for 
(3 = 1.80 on a 8^ X 4 lattice. The solid vertical lines represent the first two 
RW critical lines {fij = tt/8 and fij = Sn/S). 
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Figure 4: As in Fig. 3 for /? = 1.35. 
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Figure 5: Polyakov loop (above) and chiral condensate (below) vs fi for 
j3 = 1.90 on a 8^ X 4 lattice. Numerical data for imaginary chemical potential 
(squares) are fitted by a fourth order polynomial (long-dashed line). The 
"rotated" polynomial, whose upper and lower bounds are given by the solid 
lines, is compared with data obtained for real chemical potential (diamonds). 
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Figure 7: As in Fig. 5 for /3 = 0.90. 
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